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ABSTRACT 

We study hadronic models of broad-band emission of jets in radio-loud active galactic nuclei, 
and their implications for the accretion in those sources. We show that the models that account 
for broad-band spectra of blazars emitting in the GeV range in the sample of Bottcher et 
al. have highly super-Eddington jet powers. Furthermore, the ratio of the jet power to the 
radiative luminosity of the accretion disc is ~ 3000 on average and can be as high as ~ 10 5 . 
We then show that the measurements of the radio core shift for the sample imply low magnetic 
fluxes threading the black hole, which rules out the Blandford-Znajek mechanism to produce 
powerful jets. These results require that the accretion rate necessary to power the modelled 
jets is extremely high, and the average radiative accretion efficiency is ~ 4 x 10~ 5 . Thus, if 
the hadronic model is correct, the currently prevailing picture of accretion in AGNs needs 
to be significantly revised. Also, the obtained accretion mode cannot be dominant during the 
lifetimes of the sources, as the modelled very high accretion rates would result in too rapid 
growth of the central supermassive black holes. Finally, the extreme jet powers in the hadronic 
model are in conflict with the estimates of the jet power by other methods. 

Key words: acceleration of particles-galaxies: active-galaxies: jets-ISM: jets and outflows- 
quasars: general-radiation mechanisms: non-thermal. 


1 INTRODUCTION 


The most commonly considered models of broad-band emission 
of radio-loud jets in AGNs utilize emission of relativistic electrons 
and positrons via synchrotron and Compton processes. These lep- 
tonic models have been highly successful in reproducing the spec¬ 
tra and time variability of blazars and radio galaxies, though there 
are some phenomena, such as their extremely short variability time 
scales_Gn some cases down to a few minutes ( lAharonian et alj2007l : 
I Albert et al. 2007), that are not readily explained by them. 


The alternative, and quite popular, model is_based on hadronic 
(or lepto-hadronic) processes (e.g., Mannheim & Biermannl 11993: 


[Miicke & Protheroett200ll : Aharonianl l2000i ). It utilizes emission of 

extremely relativistic ions, mainly due to proton-synchrotron and 
photo-pion production processes, to explain the high-energy emis¬ 
sion from jet-dominated AGNs. The latter process gives rise to cas¬ 
cades of electrons and positrons, also emitting via the synchrotron 
and Compton processes. The peak of the high-energy emission of 
the broadband spectral energy distribution of blazars, at ~ 0.1-10 
GeV, is then due to these processes, mostly proton synchrotron. 


We stress that the name ’hadronic’ refers only to the emission 
processes, and the ’leptonic’ jets in AGN models also contain ions, 
though with relatively low energies, preventing them from radiating 


efficiently. Still, those hadrons in the leptonic models usually dom¬ 
inate ffiekineticluminosity o f jets (e. g., Ghisellini et al.ll2014l) . 

ISikora et all (2009 ) and ISikoral 1201 ll) have shown that the 
hadronic processes are very inefficient. This implies that, if the to¬ 
tal jet power is approximately limited by the Eddington luminosity, 
the hadronic model can be ruled out in many cases. On the other 
hand, iBottcher et all 120131 ) (hereafter B13) have successfully ap¬ 
plied hadronic models to a sample of radio-loud AGNs circumvent¬ 
ing this constraint by allowing the jet power to be highly super- 
Eddington. Here we discuss consequences of this supposition. 


2 ANALYSIS OF THE SAMPLE OF B13 

We study here the sample of 12 blazars of B13. Their broad-band 
spectra were fitted by them by leptonic and hadronic models, and 
we consider here the jet powers obtained by B13 for the latter. B13 
provide the redshifts, z, the apparent superluminal velocity, /l app , 
and the accretion luminosity, L a „„ , see Table 1. Six of those objects 
were also present in the sample of Zamanin asab et al] d20 1 4) (here¬ 
after ZCST14), for which those authors provide estimates of the 
black-hole mass, M,/? app , and L acc based on literature. The values 
of yS app agree well between ZCST14 and B13, while the values of 
L acc agree to within a factor of ~2, which is a satisfactory agree- 
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Figure 1. The ratio of the total jet power to (a) the Eddington luminosity 
and (b) the accretion luminosity as functions of the Eddington ratio. The 
dashed lines show the geometric averages. 


ment given the necessary approximate character of those estimates. 
However, the exception is W Comae, for which the L acc estimate 
used by B13 is 40 times higher than that given in ZCST14, which 
is the upper limit o f [Ghisellini et alj OOld) . hereafter G10. B13 
used the estimate of lXie et alj|2008l) , who refer in turn to two other 
papers, which do not seem to give that information. Thus^we use 
here the value of G10. B13 quote an estimate of L acc of lxie et~akl 
d2008h for OJ 287, which is 7 times higher than the upper limit of 
G10; here we use the latter. For two objects (S5 0716+714 and 3C 
66A), B13 give no estimates of L acc , and we use the upper limits 
given by G10. In the case of 3C 66A, there is also an uncertainty 
of its redshift, and B13 used an estimated redshift of z = 0.3 based 
on EBL-corrected SED modelling iAbdo et al.l201ll ). which is also 
marginally consistent with the spectro scop ically determined lower 
limit on the redshift bv lFumiss et al.lj2013l) . Since we use the esti¬ 
mates of the jet power from B13, we adopt their redshift value and 
rescale the upper limit of L acc of G10, who used z = 0.444, by the 
ratio of D 2 L , where D L is the luminosity distance. We assume the 
same cosmological parameters as B13, Sl\ = 0.7, fi ra = 0.3 and 
Ho = 70 km/(s Mpc). 

For the objects not in the sample of ZCST14, we find mass es¬ 
timates from the literature, in particular from G10. In some cases, 
there is a considerable uncerta inty. For OJ 287, we use 6.2x 10 8 Mq 
from 1 Wang. Luo & Hoi ( 2004}) . However, if the system is indeed a 
binary black-hole system, the mass of the less and more massive 
component is estimate d as 1. 4 x 10 8 Mf?) a nd 1.8 x 10 10 Mq, re¬ 
spectively l lValtonen. Ciprini & Lehtoll2012l) . Accretion would be 
dominated by the more massive black hole. The adopted values of 
L acc and M for all sources studied here are given in Table 1. 

B13 give the j et poweif] estimated using the energy content and 
for one jet only, see their equations (3-5 ). However, the jet power, 
Pj, is the enthalpy flux (e.g.. lLevinsorj|~2006l) . and the counterjet 
should be included in the energy budget. Thus, we multiply their 
values by 8/3, assuming the protons are relativistic (and thus their 


1 The entry for Lp (the power in protons) of OJ 287 in B13 is a typo, it 
should be 0.083 rather than 8.3. 
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Figure 2. The dimensionless magnetic flux 0 bh as a function of the Ed¬ 
dington ratio. The dashed line show the geometric average of 0 bh - 1-9. 
The model of efficient jet f ormation via black- hole spin energy extraction 
JBlandfor^&Znaje^ 19771) predict 0bh ~ 50 jTchekhovskov et al]l201 ll : 
Mc Kinney et all2012 ), which is much higher than the obtained values. 
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equation of state can be described by an adiabatic index of 4/3), 
see Table 1. We neglect a possible contribution of cold protons, and 
consider only those in the power-law distribution assumed to give 
rise to the observed spectrum, following the fits of B13. We then 
compare Pj to both L E (calculated for an H fraction of X = 0.7) and 
L acc in Fig. [J] We find (Pj/L E ) - 85, ( Pj/L acc > - 2660, with the 
standard deviations corresponding to factors of =H0 and ^8.4, re¬ 
spectively. (Hereafter the symbol (.) denotes a geometric average, 
i.e., average of logarithms, and the given standard deviation corre¬ 
sponds to a multiplicative factor.) Except for W Comae and OJ 287, 
the distribution of Pj/L E is relatively uniform. As a consequence, 
there is an anti-correlation between fj/L acc and L acc /L E , but all the 
obtained values are still very high, spanning Pj/L- dcc ~ 50-10 5 . Fur¬ 
thermore, some of our values of L acc are upper limits, and the cor¬ 
responding true values of Pj/L dCC can be even higher. 

In general, we can write Pj = ejMc 2 , where M is 
the accretion rate and e ) S; 1.5, the ratio of the jet 
power to the available accretion power, which we will call 
here the jet formation efficiency. The approximate maxi¬ 
mum corres ponds ^Tchekhovskoj^jjara Yan & McKinney! 1201 ll : 
McKi nney, Tchekhovskov & Bland(orJ|20J2|) to efficient ext rac- 
tion of the black-hole spin energy iBlandford & Znaiekll 19771) . in 
which case the jet power may exceed Me 2 . However, attaining 
such a high jet formation efficiency requires the accretion flow 
to possess a strong poloidal magnetic field (Tchekhovsko v et al.l 
l201ll : [McKinney et al. 201 j), namely to form a magnetically ar- 


rested d isc (MAD, Narayan, Igumenshche v & Abramowic j|2003l : 

see also iBisnovatvi-Kogan & Ruzmaikinlll976l) . The required en¬ 
ergy density corresponds then to the magnetic flux threading a 
rapidly spinning black hole of ®bh = 0bh(M c) l/2 r g , with bh — 50, 
which corresponds to the magnetic field strength at the horizon of 
B 2 /8n ~ 100/m?ZpC 2 /(<r T r g ), where m = Me 2 /L E , r % = GM/c 2 is 
the gravitational radius, m p is the proton mass, and cr T is the Thom¬ 
son cross section. Since the magnetic flux is conserved in the jet, 
it can be measured at large scales, e.g. , through the radio-core shift 
effect jLobanovll 19981 : lHirotanj||2005l) . Values of 0 bh ~ 50 were 
found for a large sample of blazars and radio galaxies by ZCST14 
assuming Me 2 = L acc /e r for their adopted accretion radiative effi¬ 
ciency of e r = 0.4. However, given that Pj » L acc in the hadronic 
model of B13, we have to estimate M instead from the jet power, 
strongly dominating the energy budget, i.e.. Me 2 — PjAj » L acc /e r . 

To calculate the values of 0 bh implied by the hadronic model, 
we use the radio core shift, Ad, measurements_between 8 and 15 
GHz for the sources in our sample from IPushkarev et alj J2012h. 
which we give in Table 1. We apply a correction of Zdziarski et al.l 
120151) to the power of the (1 + z) factor in the formula for the mag- 
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Table 1. The main parameters of the sample. FSRQ - fiat spectrum radio quasar. LBL, IBL and HBL refer to blazars which synchrotron spectrum has the peak 
at low, intermediate and high frequencies, respectively; v and vF v give the estimated position of the peak of the GeV spectrum approximated as a power law, 
and Pj is the total jet+counterjet power in the hadronic model ofB13.SeeB13 for the remaining parameters of the sources. 


Object 

type 

Z 

Lacc 

M 

Lace / Le 

A0 

V 

vF v 

Tmin 

Tmax 

P> 




erg s _1 

Mq 


mas 

10 22 Hz 

10~ 10 ergcm _2 s -1 


10 s 

10 48 erg s -1 

0219+428 (3C 66A) 

IBL 

0.3 

4.80X10 43 

4x 10 s 

8.ixi0' 4 

0.073 

10 

1 

1 

1.2 

3.2 

0235+164 (AO) 

LBL 

0.94 

3.40X10 44 

4x10 s 

0.0057 

0.006 

30 

7 

1 

4.3 

27 

0420-014 (PKS) 

FSRQ 

0.914 

3.02X10 46 

2.57x10 s 

0.97 

0.256 

5 

1 

10 3 

0.43 

1.4 

0528+134 (PKS) 

FSRQ 

2.07 

1.70X10 47 

1.07x10 s 

1.07 

0.150 

5 

2 

1 

1.1 

117 

0716+714 (S5) 

LBL 

0.31 

1.80X10 44 

lx 10 s 

0.012 

0.125 

70 

1 

1 

2.7 

0.69 

0851+202 (OJ 287) 

LBL 

0.306 

1.50X10 44 

6.2x10 s 

0.0016 

0.051 

5 

0.5 

1 

1.0 

0.23 

1219+285 (W Comae) 

IBL 

0.102 

1.50X10 43 

5x10 s 

2.0x 10~ 4 

0.047 

80 

0.5 

1 

1.9 

0.059 

1226-023 (3C 273) 

FSRQ 

0.158 

1.30X10 47 

6.59x10 s 

0.13 

0.017 

2 

10 

10 3 

0.43 

67 

1253-055 (3C 279) 

FSRQ 

0.536 

2.00x 10 45 

8x10 s 

0.017 

0.051 

10 

1 

10 3 

0.63 

9.3 

1510-089 (PKS) 

FSRQ 

0.36 

1.12X10 46 

1.58x10 s 

0.48 

0.151 

8 

5 

10 3 

1.1 

6.7 

2200+420 (BL Lac) 

LBL 

0.069 

1.51X10 45 

1.70x10 s 

0.060 

0.052 

10 

0.3 

1 

1.9 

26 

2251+058 (3C454.3) 

FSRQ 

0.859 

7.24x 10 46 

4.90x10 s 

1.00 

0.159 

10 

10 

10 3 

1.1 

96 


netic field strength at 1 pc, B\, and then use the formula for the jet 
magnetic flux, Oj, of equation (5) of ZCST14. From Oj = fl> BH we 
obtain 0 bh. which we plot in Fig. [2] for e } = I and the dimension¬ 
less black-hole angular momentum of 1. We find (0 bh) —1.9 with 
the standard deviation of a (multiplicative) factor o f ^ 5. These re¬ 
sult can be compared to another formula for Oj of lZdziarski et alj 
I20151) . which takes into account effects due to the small jet opening 
angles in the radio-emitting region (and the consequent low value 
of the magnetization parameter), and due to transverse averaging 
of the magnetic field. That formula yields values lower by approxi¬ 
mately V2 from those above. Thus, our results appear to give robust 
approximate estimates of the magnetic flux (under the assumption 
of the hadronic jet model). Given this and that all objects in the 
sample h ave </i R h 50, we can rule out the efficient version of the 
model of lBlandford & Znaiekl d 1977h for the sample. (We note that 
this model is likely to work for the studied sources if the leptonic 
model is adopted, as in ZCST14.)_ 

In other jet formation models (e.g. jBlandford & PavneilT982l : 
ICoughlin & Begelmanl 120141) . the jet formation efficiency has to 
be < 1, and it is usually ej <K 1. Thus, the accretion pro- 
c ess is even more high ly su per-Eddington. For such accretion, 
ICoughlin & Begelmanl 1 12014b make a rough estimate of ej = 0.1, 
which we also adopt. Then, (Mc 2 /L E ) ~ 10 3 , i.e., the accretion is 
indeed highly supercritical. The radiative efficiency of the accretion 
disc in this case is very small, (L acc /Mc 2 ) =; 4 x 10 -5 . 

After the calculat ions presented in this paper had been com¬ 
pleted, ICerruti et al.l < 120151) presented an analysis of 5 high- 
frequency peaked BL Lac objects (HBLs), whose y-ray peak is lo¬ 
cated around ~ 1 TeV rather than at ~ 1 Ge V as typic al for LBLs and 
IBLs considered in B13 and in this work. lCerruti et all d2QI5l) found 
that it is possible to find hadronic-model fits with sub-Eddington jet 
powers for those particular objects. A detailed analysis of those ob¬ 
jects is outside th e scope of th is Letter, but we briefly address a 
few issues. ICerruti et alj d2015l) took into account secondary prod¬ 
ucts of hadronic processes in more detail than B13. However, they 
assumed steady state primary proton and electron distributions as 
fixed power laws and followed the evolution of particle spectra, 
using the respective kinetic equations, only for the secondary par¬ 
ticles, while B13 used the electron and proton kinetic equations for 
all particles, which is the self-consistent approach. It is not clear to 
us how this might affect the jet powers. Also, the jet powers given 
in ICerruti et al. 1 120 151) are for one jet only, neglect the pressure con¬ 


tribution, and assume the jet Lorentz factor, T, to be one-half of the 
Doppler factor, S , while T = S for the inclination angle i = 1 /T, as 
favoured on statistical grounds. Thus, for comparison with the re¬ 
sults here, their jet powers need to be multiplied by a factor of 32/3. 
Still, we do not exclude that a sub-Eddington jet powers can be ob¬ 
tained for some objects, as it is, e.g., the case for W Comae in our 
sample. Generally, due to the substantially lower total luminosities 
and at least equally large (if not larger) black-hole masses of HBLs 
compared to low-frequency peaked blazars (especially FSRQs), it 
seems less problematic to produce hadronic model fits with sub- 
Eddington jet powers for HBLs. 


3 THE MINIMUM JET POWER IN THE 
PROTON-SYNCHROTRON MODEL 

The above results are for the model fits of B13 which, due to the 
large number of parameters, may generally not be the only possible 
hadronic model representation of the chosen SEDs. We may there¬ 
fore ask whether another set of fits of the hadronic model could 
have lower total power. We can answer this question by finding 
the minimum possible power for a given proton synchrotron flux, 
found to dominate the overall model spectra of the objects in the 
sample, see fig. 9 of B13. We use the method of lZdziars 3B 
(hereafter Z14), who calculated the minimum jet power for a given 
electron synchrotron flux. We can adapt those calculations to the 
proton synchrotron case by replacing the electron mass, m t , by 
m p in the relevant formulae of Z14 (including physical constants). 
Equivalently, <r T , the critical magnetic field, B cr , the dimensionless 
photon energy in the jet frame, P, and the jet-frame flux, L' € , (in 
the notation of Z14) need to be multiplied by (m p /m e L 2 , (m p /«7 e ) 2 , 
(/7z p //7t e ) 1 and m p /m e , respectively. Correspondingly, the constant 
a e of equation (33) of Z14 needs to be multiplied by (m p /m e ) 3/2 , and 
the power in relativistic particles, P c , and in the rest mass, Pj, need 
to be multiplied by (m p /m e ) 5/2 and m p /m e , respectively. The mini¬ 
mum total jet power and the corresponding magnetic field strength 
of equation (36) in Z14 need to be multiplied by (m p /m c ) l0/7 and 
(m p /m e ) 5/7 , respectively. 

The method of Z14 assumes a local synchrotron spectrum 
which can be well represented by a power law in some range of 
frequencies, and it uses one monochromatic measurement of the 
flux. Here, we use an estimate of the flux around the maximum of 
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Figure 3. (a) The ratio of the total jet power estimated using the magnetic 
field strengths of B13 and the method of Z14 to the corresponding power 
given in B13, as a function of the Eddington ratio. The dashed line corre¬ 
sponds to equal powers, (b) The ratio of the minimum power to that calcu¬ 
lated using the method of Z14 with the magnetic field strengths of B13. The 
dashed line corresponds to the geometric average. 


the vF v spectrum, using spectral plots in fig. 9 of B13. We list our 
adopted values in Table 1. This implies that the parameter e max of 
Z14 is = 1. However, we modify the method in order to take into 
account that y m i„ is specified in the model of B13. In the case of ion 
acceleration, one would naturally expect y m i n ~ 1, unlike the case 
of electrons, which can be preheated to (y min ) » 1 (with acceler¬ 
ation only proceeding out of this distribution). However, B13 used 
ymin = 10 3 for some objects, see Table 1. To account for this, we 
define e min = ByPj{B a m 2 e') (cf. Z14). Since B is itself calculated 
by the minimization, this requires a simple iterative calculation. 

We first check how well the above method reproduces the jet 
powers given in B13. For this, we use the values of B given in B13 
instead of those corresponding to the minimum power. We show 
the results in Fig. (3ja). Some disagreement is expected given that 
in the present estimate we take into account only the proton syn¬ 
chrotron emission and neglect cooling. The adiabatic cooling is, in 
fact, dominant in AO 0235+164 and W Comae, in which case the 
proton-synchrotron estimate gives values ~200 and 50 times too 
low, respectively. Apart from this, the overall agreement is good, 
within a factor of 2 for 7 for the remaining objects, see Fig.[3ja). 

We then apply the minimization method to the sample of B13. 
We have found the minimum powers are on average lower by a 
factor of few compared to the corresponding estimates above, with 
the geometric average being ^ 0.23, see Fig. [3])b). This reduction 
is due to the minimum corresponding to an equipartition between 
the power in the relativistic ions, P,, and in the magnetic field, P B , 
which may not be able to produce a fit to the entire SED. Indeed, 
B13 found P[ » P B in most of their fits. This departure from 
equipartition naturally results in higher total jet powers. Still, the 
overall possible reduction of the jet power due to changing the fit 
parameters is by at most a factor of a few. 


4 ASTROPHYSICAL IMPLICATIONS 


In Section [2] we found that the jet powers of the AGN sample 
of B13 exceed their radiative luminosities by very large factors, 
with Pj/Lacc ~ 50-10 5 and (Pj/L acc ) ~ 3000. Then, we showed in 
Section [3] that the values of the jet power obtained by B13 in the 
hadronic model cannot be significantly reduced by changing the 
fit parameters, as they are, on average, within a factor of a few of 
the minimum jet powers estimated using the observed y-ray spec¬ 
tra for the radiative mechanism dominant in this model, namely 
proton synchrotron. These findings imply that the jet formation 
efficiency greatly exceeds that found in leptonic models, which 
is already quite high (see, e.g.. iGhisellini et all l2014h . We have 
also found that the most efficient jet formation mechanism known, 
based on the Blandford-Znajek mecha nism extracting black-hole 
spin with a magnetica lly arrested disc iTchekhovskov et al]|201ll : 
McKi nney et alj2012l) . is ruled out by the magnetic fields measured 
at the pc distance scale. Since the jet power has to be then derived 
entirely from accretion, the implied accretion rates are huge and the 
radiative efficiencies of the accretion disc must be tiny. 

At the inferred accretion rates, (Mc 2 /Le) ~ 10 3 , accre¬ 
tion is supercritical (ruling out optically-thin radiatively inefficient 
models, which can co rrespond only to Me 2 /Le <K 1, see, e.g., 
lYuan & Naravanll2014l) . The angle-averaged accretion-disc lumi¬ 
nosity in this case is small, corresponding to L ~ L E and a ra - 
diative efficiency is ~ 10 3 dSikoralll98ll : ISadowski et al1 l20T4h . 
However, most of that luminosity emerges through axial funnels, 
where the observed flux is super-Eddingto n, corresponding to L ~ 
10L e JSikorall 198 ll : ISadowski et alj|2014h . Our sample consists of 
blazars seen close to on-axis, on average at i ~ 1 /Fj, and the 
accretion-disc-emission funnels have opening angles greater than 
that ( ISadowski et alJl2014t) . Thus, we would expect to see super- 
Eddington accretion-disc luminosities, while, in fact, they are all 
;$ L E , with 8 out of 12 objects having L. iCC S 0.1L E . Thus, the 
hadronic emission model for the jets is inconsistent with the stan¬ 
dard accretion theory. 


The inferred extreme accretion rates present also a major 
problem in light of results of studies of supermassive-black-hole 
growth. The e-folding time by which the black-mass would in¬ 
crease due to accretion, M/M, is = 4 x 10 5 yr at {Me 1 /Le) ^ 10 3 
found in Section[2] This is a very short time compared to estimated 
lifetimes of active phases of radio_sources, e.g., = 2 x 10 s yr for FR 
IIs dAntognini. Bird & Martinill2012i and references therein). Cor¬ 
respondingly, the average radiative accretion efficiency we found 
(under the assumptions of the hadronic model) of (L acc /Me 2 ) ^ 4x 


10 5 i s <K than the average accretion efficiency of ~ 0 .1-0.3 ij SoItan 


1982; iMarconi et aid 120041 : Isilverman et aid 120081 : ISchulze et al 


2013). Thus, blazars radiating via the hadronic emission mecha¬ 


nism have to represent at most a small fraction of accreting super- 
massive black holes and/or such extreme accretion episodes must 
be extremely short-lived, representing only a duty-cycle of order 
~ 10~ 4 (which is in conflict with the fact that the studied objects 
have quite average properties). 


Finally, the obtained extreme powers are in conflict 
with studies o f the jet power base d on radio lobes and 
X-ray cavities dM erloni & Heinz! 12 0071: ICavagnolo et ahl I201C ; 


Nemmen et a 1 J 20 T^T Russell et alJ2013l : lGodfrev & Shabalal2013 : 

Shabala & Godfrey! 120131) . Those studies indicate jet powers at 
most moderately exceeding the Eddington luminosity, and even in 
the most luminous sources they exceed the accretion-disc luminos¬ 
ity only by a factor of ^ 10. On the other hand, those powers are in 
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an overall agreement with the powers estimated in leptonic models, 
see, e.g. jGhisellini et ahl | |2014|) . 

Summarizing, we have found major difficulties in reconciling 
the jet power requirements of hadronic blazar models with (a) ob¬ 
served accretion-disc luminosities, (b) accretion rates inferred from 
supermassive-black-hole growth, and (c) jet powers inferred from 
radio lobes and X-ray cavities. While the problems (b) and (c) can 
be, in principle, circumvented, if we apply the hadronic model to a 
small number of sources and assuming a short duty cycle, the prob¬ 
lem (a) requires a change of the accretion paradigm for the very 
sources fitted by the hadronic model. We conclude that our results 
represent a strong argument against the applicability of hadronic 
models to blazars. 


5 SUMMARY 

We have presented an analysis of the jet powers required by 
hadronic model fits to the SEDs of blazars, based both on detailed 
numerical modelling by B13, and on minimal jet power require¬ 
ments in a proton-synchrotron interpretation of the y-ray emission, 
adopting the methodology of Z14. Our main results are as follows. 

We show that hadronic models of B13 that can account for 
broad-band spectra of blazars emitting in the GeV range, have jet 
powers of P- 3 ~ 10 2 Le, approximately independent of their accre¬ 
tion Eddington ratio, which spans L dcc /L E - 10 3 —1 for the studied 
sources. The ratio of the jet power to the radiative luminosity of the 
accretion disc can be as high as ~ 10 5 for the studied sample, and 
(Pj/Ltcc) ~ 3000. 

Furthermore, we show that the available measurements of the 
radio core shift for the sample imply low magnetic fluxes threading 
the black hole, which rules out the Blandford-Znajek mechanism 
with efficient production of jets for hadronic blazar models. 

These results require that the accretion rate necessary to power 
the modelled jets is very high, compared to the accretion radiative 
output. The average radiative accretion efficiency is (L dCC /Mc 2 ) - 
4 x 10~ 5 for the studied sample. 

A major problem for the the hadronic model is presented by 
their required highly super-Eddington jet powers. This, in turn re¬ 
quires highly super-Eddington accretion rates, at which the ob¬ 
served luminosities would be much higher than those actually seen. 
If the model is correct, the currently prevailing picture of accretion 
in AGNs needs to be revised. 

The extreme jet powers obtained by B13 are in conflict with 
the estimates of the jet power by other methods. The obtained 
accretion mode cannot be dominant during the lifetimes of the 
sources, as the modelled very high accretion rates would result in 
much too rapid growth of the central supermassive black holes. 
Thus, if applicable, this accretion mode can only be present with 
a very small duty-cycle in the black-hole evolution, and thus can be 
applied only to a very small fraction of the radio loud sources. 
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